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ABSTRACT 

As the brightest 0-type X-ray source in the Carina nebula, HD 93250 (04 Ill(fc)) is X-ray 
overluminous for its spectral type and has an unusually hard X-ray spectrum. Two different 
scenarios have been invoked to explain its X-ray properties: wind-wind interaction and magnetic 
wind confinement. Yet, HD 93250 shows absolutely constant radial velocities over time scales 
of years suggesting either a single star, a binary system seen pole-one or a very long period 
and/or highly eccentric system. Using the ESO Very Large Telescope Interferometer, we resolved 
HD 93250 as a close pair with similar components. We measured a near-infrared flux ratio of 
0.8 ± 0.1 and a separation of 1.5 ± 0.2 x 10~^ arcsec. At the distance of Carina, this corresponds 
to a projected physical distance of 3.5 A.U. While a quantitative investigation would require 
a full characterization of the orbit, the binary nature of HD 93250 allows us to qualitatively 
explain both its X-ray flux and hardness and its non-thermal radio emission in the framework 
of a colliding wind scenario. We also discuss various observational biases. We show that, due to 
line-blending of two similar spectral components, HD 93250 could have a period as short as one 
to several years despite the lack of measurable ipdial velocity variations. 



Subject headings: Stars: early-type - Stars: individual: HD93250 - Stars: massive - X-rays: individual: 
HD93250 - Techniques: interferometric - Radiation mechanisms: non-thermal 



1. Introduction 

Massive stars of spectral-type O are typically 
soft and bright X-ray sources. Their X-ray emis- 
sion is usually attributed to shocks in the dense 
region of the winds that grow from the intrinsic 
instabilities of the line driving mechanism. In the 
presence of a relatively strong magnetic field, the 
structure of the wind might be controlled by mag- 
netic force lines and, at least in case of a bipo- 
lar field, the wind might be confined towards the 
equator, producing a typically harder X-ray emis- 
sion compared to that expected from embedded 
wind shocks. In binary systems, an additional 
X-ray emission can also be generated by the su- 
personic collision of the two winds. One of the 
most comprehensive surveys of the X-ray proper- 
ties of massive stars has recently been completed 
in the framework of the Chand ra Carina Complex 
Project (Townslev et al.ll2011l ). that has provided 
X-ray spectroscopy of 68 0-type stars in the Ca- 
rina nebula. Yet not all sources are equally-well 
understood. 

Located at 7.5' NNW from Eta Carinae in 
the Trumpler 16 com plex, HD 93250 (04 Ill(fc); 



where the collision between the wind of the two 



Walborn et all I2OIOI ) is the brightest X-ray O 



star in Carina and has puzzled astronomers for 
over two decades. Assuming that it belongs to 
the Carina association, HD 93250 is too bright 
for its spectral type. Its X-ray lumino sity, 



1.5 X lO^^ergs-i (jCagne et all l201l[ ). 



is 



a factor two to three larger than expected from 
the canonical Lyj — Lhoi relation for O-type stars 



([Sana et al.l 120061 : iNaze et all 120111 ). Its X-ray 



emission is variabl e at the ^20 - 50% l e vel on a time 
scale of months ( Rauw et al. 2009t ICagne et al 



201 1[ ) and its X- ray spectrum is ex tremely hard 
(kT ~ 2.3 keV, lOagne et al.ll2nilh . HD 93250 
is also a known non-thermal (NT) radio emitter 



( Leitherer et al.lll995r ) 



Most of the puzzling properties of HD 93250 
could be explained as part of a binary scenario 
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stars produces both an X-ray excess ()Luo et al 



L1990I) and an NT radio emission component in the 
form of synchrotron radiation in the presence of 
a po pulation of relativistic electrons (|De Becker 



20071 ). The latter point is especially important in 



the sense that it provides the opportunity to in- 
vestigate particle acceleration processes in massive 
star environments, in addition to evolved objects 
such as supernova remnants or binaries containing 
a compact object. For the wind collision scenario 
to produce enough X-ray and synchrotron radio 
flux, the winds from both components need to be 
of similar strength, and the companion is thus ex- 
pected to be a bright O star. Yet, both high an - 
gular resolution observat ions (iNelan et al" 2004) 
and RV campaigns (e.g., iRauw et al.ll2009l ) have 
failed to reveal a binary nature of the system. This 
has raised the question whether massive single 



stars could produce NT emission (Ivan Loo et al 



20061) . As an alternative scenario, iRauw et al 



( 2009t ) suggested that the X-ray over-luminosity 
cou ld originate from a magnetically confined wind 
but lCagne et al.l (|201ll) noted that the strong non- 
thermal emission argued against such a scenario 
and called for further observations. 

As discussed by ISana fc Evand ( 2010l ). long 
base-line interferometry is ideally suited to iden- 
tify binaries in a separation regime that cannot 
be reached by either spectroscopy, typically lim- 
ited to semi- major axes up to a few A.U., nor by 
other high resolution techniques, typically limited 
to a minimum angular separation of 20 mas, corre- 
sponding to ~40 A.U. at the distance of the Carina 
complex. 

In this context, we observed HD 93250 using the 
Very Large Telescope Interferometer (VLTI) in an 
attempt to prove the binary nature of the system 
(Sect. [2|) and to perform a first characterization of 
the physical size of the orbit (Sect. [3]). 

2. Interferometric measurements 
2.1. Observations and data reduction 



We observed 



(jHaguenauer et al 



HD 93250 
2OIOI) 



at the ESO/VLTI 

with the Astronomical 

Mult i-BEam combineR (AMBER, IPetrov et aTl 
20071 ) on December 27 2010. Observations have 
been made in configuration AO-KO-Gl, providing 
projected baselines ranging approximately from 



2 



50 m to 120 m. We recorded spectrally dispersed 
fringes between 1.55 — 2.4 at a spectral resolv- 
ing power (R) of 40. 

At VLTI the fast telescope guiding (STRAP) is 
performed in the V-band. For stars fainter than 
y 11 mag, the quality of the guiding is de- 
graded and the amount of flux injected into the 
instrument is dramatically reduced. The offered 
limiting magnitude of AMBER {H — 5.5 mag) 
has been measured on late- type stars, for which 
H ^ 6 mag typically corresponds to F 11 mag. 
The offered magnitude thus results from the com- 
bination of the low flux limit in the observation 
band with the poor performance of the guiding 
for stars faint in V . 

Early-type stars are intrinsically 0.7-0.8 mag 
bright er in the V^-band com pared to in the H- 
band ( Martins fc Pled I |2006[ ) although the appar- 
ent V — H color depends on the extinction prop- 
erties. For stars as bright as HD 93250 in the 
y-band (V = 7.5 mag) the telescope guiding per- 
forms very well. The fraction of flux injected into 
AMBER is optimal and the sensitivity of the ob- 
servations is only limited by the intrinsic flux of 
the object in the observation band. As a conse- 
quence, HD 93250 [H — 6.7) could be observed 
with the 1.8m Auxiliary Telescopes (ATs) despite 
the fact it is significantly fainter than the offered 
limiting magnitude. 

Calibrated visibility and closure phase values 
were computed using the la test public version 
of the amdlib pa ckage fv3.1. iTatulH et all 12007 : 
Chelli et al. 20091 ) and the yorick interface pro 



vided by the Jean-Marie Mariotti Centeifl]. The 
main parameters of this pipeline are the criteria 
and the threshold used to reject individual frames. 
Based on our experience, we kept the 50% best 
frames sorted by signal-to-noise ratio. We also dis- 
carded the fringes recorded with an optical path 
difference larger than 20 /zm. For each spectral 
channel, we estimated the transfer function for 
square visibilities and closure-phase by averaging 
the measurements done on the nearest calibration 
star. This transfer function was divided and re- 
moved from the visibilities and closure-phases ob- 
tained for the scientific target. The dispersion 
between the consecutive points obtained for each 
observation of a calibration star is an estimate 



of the transfer function instability. This insta- 
bility was quadratically added to the uncertainty 
on the transfer function. It is the dominant un- 
certainty term for visibilities, but is negligible for 
closure phase. The final uncertainty on the trans- 
fer function was propagated to the uncertainties 
on the calibrated quantities obtained for scien- 
tific targets. Note that these procedures are now 
fully embedded in the amdlib package. Concern- 
ing the wavelength calibration, we estimate the 
accuracy of our wavelength table to be 0.05 /zm 
(that is about 4%) by cross-correlating our ob- 
served s pectra with an atm ospheric model as de- 
tailed in lKraus et al.l ( 20091 ). 



2.2. Modeling 

The absolutely calibrated visibilities and clo- 
sure phases are reported in Fig. [T] The target is 
well resolved in one baseline (KO-AO) and signifi- 
cantly less resolved on the two others (Gl-AO and 
KO-Gl). The closure phase shows a clear depar- 
ture from (or 180) degrees. As expected, it is 
impossible to fit the data with a single-disk model 
as this yields an unrealistic diameter of 3.1 mas, 
a final x^^^j ss 10 and systematic discrepancies in 
the fit. Altogether, the binary nature of HD 93250 
is conclusively revealed by these observations. 

We fitted the complete set of data (visibilities 
and closure phases) with a model composed of two 
unresolved stars with constant flux ratio over the 
H and the K bands, a valid assumption for 0-1-0 
binaries. The separation, the orientation and the 
flux ratio were the only free parameters of the fit. 
We then explored the full parameter space to find 
the best value and we obtained a very good 
fit with Xrcd ~ 0-9 J as shown in Fig. [TJ Thanks 
to the combinations of three baselines of different 
orientations and of the H and the K bands, the 
position angle of the binary is unambiguous. By 
chance, an additional "locking point" for the bi- 
nary model is provided by the visibility minimum 
observed in baseline KO-AO. 

The best fit results are summarized in Table [T] 
The exact shape of the astrometric error ellipse 
can be seen in the X^^^ plane of Fig.[T] In Table[TJ 
p is the projected separation; 9, the position angle 
measured East from North; and i2,K/-^i,K, the 
if-band flux ratio. 



-^JMMC software: http://www.jiimic.fr 
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Fig. 1. — Square visibilities (top panels) and closure phase (bottom left panel) from the AMBER observations 
of HD 93250 compared with the best binary model. The bottom right panel provides the Xrcd ni^-P with 
the flux separation of the best binary model. The best value is x^cd ~ ^.9. The red contour corresponds to 
= 4.5. 
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2.3. Additional observations 

HD 93250 has been observed two more times 
with VLTI/AMBER. On March 24 2011, we used 
the 8.2-ni telescopes (UTs) in configuration UTl- 
UT3-UT4 but we could not resolve the binary. 
This is in agreement with the very small an- 
gular separation between the components and 
the smaller baseline provided by the UTs. On 
March 29 2011, we used the ATs and we re- 
solved HD 93250 again. We measured a bi- 
nary separation of dE — 1.37 ± 0.25 mas and 
dN = 0.52±0.75 mas and a flux ratio of0.85±0.10, 
which is fully consistent with our observation of 
December 2010 reported in Table [T] Given the 
uncertainties, a rotation on the plane of the sky 
smaller than about 30° (thus 1/12 of an orbital 
cycle) would remain within the error bars. This 
second observation thus suggests a period larger 
than 3 yr, a high inclination or a very eccentric 
system. Alternatively, we might have observed 
HD 93250 at the same orbital phase (±0.5, as the 
sign of the closure phase is unconstrained in the 
observation of March 2011). 

3. Discussion 

3.1. Optical observations 

Rauw et al.l (|2009l ) measure RVs from six op- 



tical spectra of HD 93250, obtained from 1999 
to 2004, and concluded that variations are ab- 
sent, with a measur ed 1-a disper sion of 1.3 kms~^. 
We reanalyzed the iRauw et al.l data set, comple- 
mented by two additional observations obtained 
with the Fiberfed Extended Range Optical Spec- 
trograph (FEROS) and available from the ESO 
archives. The epochs of these supplementary ob- 
servations are HJD - 2,450,000 = 3365.801 (PI: 
Bouret) and 3738.828 (PI: Casassus), i.e., respec- 
tively, nine months and almost two years after the 



Table 1: Best- fit parameters of the binary model 
to the VLTI observations. 



Parameter 


Dec. 2010 


Mar. 2011 


dE (mas) 


1.44 ±0.20 


1.37 ±0.25 


dN (mas) 


0.38 ±0.45 


0.52 ±0.75 


p (mas) 


1.49 ±0.22 


1.47 ±0.35 


en 


76± 17 


69 ±27 


L2.¥i/Ll.K 


0.81 ±0.10 


0.85 ±0.10 



last I Rauw et al.l RV measurements. Despite the 
high signal-to-noise ratio (S/N) of the data (typi- 
cally about 300), we could not detect any RV vari- 
ation either. The maximum difference between the 
RVs from the different epochs is 3.6 kms~^ and 
the 1-a dispersion around the mean is 1.1 knis^^. 
We also searched for line profile variations in the 
He II and Si iv lines, with a negative result. The 
relative peak-to-peak variations of the amplitude 
and of the full width at half-maximum (FWHM) 
of e.g., the He iiA4541 line, are below 5 and 10% 
respectively, pointing again to an excellent stabil- 
ity of the HD 93250 spectrum over the years. 

3.2. Physical constraints 

HD 93250 belongs to Tr 16 in the Carin a asso- 
ciatio n. Adopting a distance of 2.35 kpc ( SmithI 
I2OO6I ). the angular separation translates to a pro- 
jected physical distance of 3.5 A.U. The second 
component appears of a quite similar brightness 
and is thus likely an O star as well. We estimate 
the mass-ratio using a mass-luminosity relation 
given by Lboi ~ M^-°^ (|Brott et al.ll201l( ). Consid- 
ering no color difference between the two compo- 
nents of the pair, thus no difference in their bolo- 
metric correction, we derive a mass ratio M2/M1 
of 0.9. 

Assuming a circular orbit and a primary mass 
Mi of 47 Mf7i as app ropriate for an 04 HI star 
(|Martins et al.ll2005a^ . a 3.5 A.U. semi-major axis 
corresponds to an orbital period of about 250 d 
and an orbital velocity of 70 kms~^. A maximum 
inclination of 2.5° is required to hide such a sig- 
nal at the 5-CT level in the noise associated the 
RV measurements. Assuming random orientation 
of binary orbits in space, one expects that such 
a low inclination is encountered with a 1:10,000 
frequency and is thus very rare. 

Alternatively, HD 93250 might be a long period, 
highly ecc entric system. W e used the method de- 
scribed in lSana et al. I (|2009h to estimate the likeli- 
hood to observe a significant RV shift as a function 
of the orbital period given the sampling provided 
by the optical observations and using no constraint 
on the orbital inclination. The probability that 
HD 93250 would have displayed a RV shift larger 
than 5 kms^^ is very close to 1.0 up to periods 
of lO'' d, then drops sharply for longer periods 
(Fig. [2]). As no such a shift is seen, this suggests 
a minimum orbital period of '-~^30 yr, which would 
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Fig. 2. — Probability that a RV shift larger than, 
respectively, 5 and 55 kms^"'^ would occur in be- 
tween any of the spectroscopic observations dis- 
cussed in Sect. 13. li as a function of the orbital pe- 
riod P. 
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Fig. 3. — Normalized composite He ii A4541 line 
profile in an 04 III-I-04 III binary with a sepa- 
ration of 50 kms~^ between the two components 
(plain lines) together with the best-fit Gaussian 
profile (dashed line) to the line profile. Gaus- 
sian noise corresponding to a S/N of 300 has been 
added to the composite profile before fitting. 



imply an eccentricity as large as 0.9 to match the 
current separation. 

There is however an additional detection bias 
that results from the blending of two equal- 
strength lines, which is not taken into account 
in [Sana et al. (2009) . Indeed the RV measure- 
ment obtained when applying single- Gaussian fit- 
ting to blended lines will remain constant (Fig. |3]) 
until evidence for a double-lined profile is seen 
(typically at ARV > 100 kms^^). To estimate 
this bias, we have generated synthetic profiles of 
04 III+04 III binaries usii ig the non-LTE a tmo- 
sphere code FASTWIND (jPuLs et all liool . for 
various RV shifts up to 200 kms~^ and for dif- 
ferent rotational broadening. Accounting for the 
typical S/N of our optical data, we have then fit- 
ted a single Gaussian profile. As expected, the 
measured RVs show no significant variations. For 
a typical rotational broadening of 100 kms~^, as 
appropriate for HD 93250, the measured ampli- 
tudes and FWHMs are changing by less that 5 
and 10% respectively for a RV difference between 
the two components of less than 55 and 70 kms~^ 
(Fig. 21). Using these new values for the detection 
thresholds, we note that the detection probability 
already starts to drop for period larger than a few 
100s days. 

As a conclusion, the fact that RV campaigns 
have failed to reveal the binary nature of HD 93250 
does not necessary imply exceptionally long pe- 
riods, high eccentricities and/or small inclination 
angles. It might also result from observational bi- 
ases affecting binaries with (quasi) identical com- 
ponents. 

3.3. Colliding winds 

Using the formahsm of Usov ( 19921 ) adapted to 
0+0 binaries, we estimate the X-ray luminosity 
that would be produced by the collision of two 
04 HI winds of equal strength to be of the or- 
der of 10^'^ ergs"^. There is enough uncertainty 
in the mass-loss rate associated with HD 93250 



(see the discussion in iGagne et al.l 120111 ) to ac- 
commodate the range of periods and separations 
discussed in the previous section. It is thus not 
possible to use the X-ray luminosity to further 
constrain the orbital properties. Assuming that 
the wind-wind collision is adiabatic, as expected 
for wide binaries, the X-ray fiux produced in the 
wind collision should scale with the separation D 
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Fig. 4. — Relative variation of the FWHM 
(/fwhm) and of the amplitude (/Ampi) of the 
He II A4541 line in an 04 III+04 III binary as 
a function of the RV separation vi — V2 between 
the two components. Three rotational broaden- 
ings are considered. 



Lx - D~} (IStevens et al.lll992l) . Given that 



Rauw et al.l (|2009l ) reported a variation of the X- 
ray flux of the order of 40%, one would thus expect 
HD 93250 to have a minimum eccentricity of 0.2. 

In a scenario where the observed radio syn- 
chrotron radiation is produced by colliding winds, 
the wind-wind interaction region should be lo- 
cated out of the bulk of the stellar wind mate- 
rial. Otherwise, the non-thermal radio photons 
suffer from significant free-free absorption and 
cannot escape the radio photosphere. Adopt- 
ing the mas s loss rate an d the te rminal veloc- 
ity given by iMartins et al. ( 2005bh (respectively 
5.62 X 10"'^ M0 yr ~^ and 3000 kms ~ ^), an d us- 



ing equations of Wright fc Barlowl (|l975l ) and 
Leitherer et"aLl (|l995l) . we estimate the radius of 



the = 1 photosphere using the same pro cedure 
as in the case of HD 168112 in ;De Beck er et al.l 
(12004). The two f requencies investigated by 
iLeitherer etall (|l995l ). i.e. 8.64 and 4.80 GHz, 



yield respectively photosphere radii of 1.4 and 
2.2 A. U. As the two components of the system are 
very similar, the wind-wind interaction region is 
expected to be located close to mid-separation, i.e. 
about 1.7-1.8 A.U. according to our results and 
assuming no significant inclination effect. Inter- 
estingly, these values are in good agreement with 
the detection of HD 93250 at 8.64 GHz, and its 



non-detection at 4.80 GHz (jleitherer et al.l[l995l ). 

We emphasize that these considerations are likely 
to be epoch-dependent. For instance, an increase 
of the stellar separation in an eccentric orbit could 
lead to the detection of the synchrotron radiation 
at 4.80 GHz as well. This scenario can only be ver- 
ified through a multi-frequency radio monitoring 
of HD 93250 on time scales covering a significant 
fraction of the orbital period. 

4. Conclusion 

Using the ESO/VLTI, we provided evidence 
that HD 93250 is an 0+0 binary system. The bi- 
nary nature gives a natural explanation of many of 
the intriguing properties of the object. In partic- 
ular, it can explain the X-ray over-luminosity and 
NT radio emission in the framework of a collid- 
ing wind scenario. The latter point lends further 
support to the so-called 'standard scenario' where 
particle acceleration occurs in the wind-wind in- 
teraction of a massive binary system. Because 
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wc sec no change in the orientation of the sys- 
tem in our two VLTI observations separated by 
three months, we conclude that HD 93250 has a 
Hkely orbital period of the order of one to several 
years. A significant eccentricity is also required to 
explain the variations observed in the X-ray flux. 

We further show that the non-detection of a 
RV signature based on eight high-resolution spec- 
tra covering a seven-year period of time does not 
require a particularly unfavorable orbital conflg- 
uration (extremely low inclination, long period 
and/or high eccentricity). It can also be ex- 
plained by observational biases resulting from line- 
blending in nearly equal-mass systems. Further 
spectroscopic and high- angular observations are 
require to characterize the orbit of the system, 
which is needed to quantitatively understand the 
X-ray and NT radio processes originating from the 
wind-wind collision. 

Together with HD 168112, HD 93250 was one 
of the only two known 0-type NT radio emitters 
for which the binary scenario was questioned. Our 
results provide further evidence that binarity and 
NT radio emission go hand-in-hand. As was the 
case for HD 93250, RV campaigns have so far failed 
to confirm the binary nature of HD 168112. Long 
base-line interferometry might thus provide a com- 
plementary diagnostic tool to test the multiplicity 
properties of HD 168112 in particular and of single 
non-thermal radio emitters in general. 

The authors warmly thank the VLTI team of 
the Paranal Observatory. CNRS is acknowledged 
for having supported this work with the alloca- 
tion of Guaranteed Time Observations. We are 
also greatful to Nolan Walborn for interesting dis- 
cussion. 
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